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ABSTRACT 



Aims. We present a study of the orbits of the triple system LHS 1070, with the aim to determine individual masses of its components. 
Methods. Sixteen new relative astrometric positions of the three components in the K band were obtained with NACO at the VLT, 
Omega CASS at the 3.5 m telescope on Calar Alto, and other high-spatial-resolution instruments. We combine them with data from 
the literature and fit orbit models to the dataset. We derive an improved fit for the orbit of LHS 1070 B and C around each other, and 
an estimate for the orbit of B and C around A. 

Results. The orbits are nearly coplanar, with a misalignment angle of less than 10°. The masses of the three components are M A = 
0.13 .. . 0.16 Mq, M B = 0.077 ±0.005 M , and M c = 0.071 ±0.004 M . Therefore, LHS 1070 C is certainly, and LHS 1070 B probably 
a brown dwarf. Comparison with theoretical isochrones shows that LHS 1070 A is either fainter or more massive than expected. One 
possible explanation would be that it is a binary. However, the close companion reported previously could not be confirmed. 
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1. Introduction 

LHS 1070 (other common names are GJ 2005, LP 881-64, 
2MASS J002444 19-2708242) is a nearby high-proper-motion 
sta r located in the sou th galactic pole region. It was observed 
bv lLeinert et al.l (1 1994b as part of their near-infrared speckle sur- 
vey for duplicity of nearby southern M dwarves. The group de- 
tected in 1993 two companions located 1.1" and 1.3" north of 
the primary by using the SHARP camera mounted at the NTT 
(La Silla, Chile) and applying Speckle imaging. In a consecu- 
tive measurement the companions exhibited an orbital motion 
around each other and around the primary while following the 
proper motion of the syste m. A fourth component D was identi- 
fied by iHenrv e t al. (1999) with the Fine Guidance Sensors on- 
board the Hubble Space Telescope (HST). A separation of only 
about 50 mas from component A was reported. Speckle interfer- 
ometric and adaptive optics measurements, however, were not 
able t o confirm LHS 1070D dLeinert et all 1200 UlSeTfahrt et al.L 
2008). This might be due to the small separation, but it is also 
possible that the detection is spurious (T. Henry, priv. comm.). 

When comparing the colours of the components, the com- 
panions appear redder than the primary ( Leinert et al.l [1994). 
First estimates already identified the two companions as stars 
close to the hydrogen burning limit that divides brown dwarves 
from main-sequence stars. This mass range is characterised by a 
strong decrease of the effective temperature and the onset of dust 
formation in the atmospheres. Photometric and spectroscopic 
measurements with the HST in the visual inde ed were rep- 
resented by model atmospheres containing dust (iLeinert et al.L 



* Based on observations collected at the European Southern 
Observatoiy, Chile, proposals number 60A-9026, 66.C-0219, 67.C- 
0354, 68.C-0539, 70.C-0476, 072.C-0022, 074.C-0637, 078.C-0386, 
380.C-0179, 382.C-0324, and 382.C-0329. 



2000). Also the low masses could be confirmed. The spectral 
classifications of the companions are M8.5V and M9-9.5V. The 
primary has an earlier spectral type of M5.5-6V. 

iReiners et al.l (|2007) found for the components B and C a 
similarily high rotational velocity of vsin/ = 16km/s, which 
is twice that of the primary. The authors conclude that the 
temperature-dependent magnetic braking was acting for 1 Gyr 
on all three components. The main difference between the oth- 
erwise very similar components B and C is the higher activity 
of component B. This might be related to the higher temperature 
or the higher magnetic flux. In the HST spectra Ho- emission 
was foun d towards LH S 1070 A and B, but not towards com- 
ponent C (ILeinert et all 12000). Recently, photometric activity o f 
the two more massive stars was identified (Almei da et alll201lb . 
Component B showed a brightness increase in the visual that is 
among the largest ever observed in a flare star. 

The very low mass companions of LHS 1070 provide the 
possibility to refine the models of (pre-)main-sequence stars 
close to the transition region to brown dwarves. It is thus very 
important to derive the exact masses of these stars. The small 
separation of the components B and C, i.e. their short orbital pe- 
riod is very well suited to allow even for this low mass stars a 
mass determination by fitting the orbital elements. The dynami- 
cal mass has the advantage of being independent from theoretical 
models. It is thus a precious probe to test evolutionary models. 

The main uncertainty for the mass determination is the dis- 
tance to the ob ject. The trigon ometric parallax of LHS 1070 de- 
rived by jvanAltenaetal] l[l995|) is 135.3 + 12.1 mas. This value 
was later refined to 129. 47 + 2.48 mas, plac ing the system at a 
distance of 7.72±0.15pc (ICosta et all l2005). The proper motion 
of LHS 1070 is 653.7 + 0.3mas/yr, corresponding to 23.9 km/s, 
and directed along the postion angle 348.3° + 0.44°. In addition, 
radial velocity measurements found that LHS 1070 is approach- 
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Table 1. Astrometric measurements of LHS 1070 B-C 



Date (UT) 



A-B 



B-C 



d [mas] 



PA[°] 



d [mas] 



PA[° 



Reference, Instrument, Proposal-No. 



1993 Jul 29 


1072 


±10 


-8.6 ±0.3 


266. 


±5. 


328.5 ±0.7 


Leinertetal. (2001) 


1994 May 1 


1085 


±4 


-7.7 ±0.1 


341. 


±5. 


346.9 ± 0.6 


Leinertetal. (2001) 


1994 Sep 15 


1092 


±9 


-7.2 ±0.4 


375. 


±4. 


352.7 ±0.61 


Leinertetal. (2001) 


1994 Sep 24 


1095 


±12 


-6.8 ± 0.42 


382. 


±5. 


353.9 ± 0.67 


Leinertetal. (2001) 


1995 Jan 17 


1094 


±6 


-6.2 ±0.3 


400. 


±8. 


356.1 ±0.4 


Leinertetal. (2001) 


1995 Jul 9 


1119 


±11 


-6.4 ±0.1 


439. 


±4. 


3.2 ±0.1 


Leinertetal. (2001) 


1995 Jul 14 


1102 


±3 


-6.1 ±0.1 


436. 


±1. 


3.7 ±0.2 


Leinertetal. (2001) 


1996 Jan 16 


1124 


±2 


-5.6 ±0.3 


459. 


±3. 


9.0 ±0.3 


Leinertetal. (2001) 


1996 Aug 22 


1161 


±5 


-5.0 ±0.1 


465. 


±3. 


14.9 ±0.1 


Leinertetal. (2001) 


1996 Sep 27 


1157 


±8 


-4.7 ±0.1 


468. 


±4. 


15.8 ±0.1 


Leinertetal. (2001) 


1997 Jul 15 


1235 


±14 


-3.5 ± 0.23 


458. 


±6. 


23.4 ± 0.6 


Leinertetal. (2001) 


1997 Aug 25 


1243 


±7 


-3.5 ± 0.23 


450. 


± 15. 


25.4 ±1.0 


Leinertetal. (2001) 


1997 Nov 17 


1223 


±5 


-2.9 ±0.1 


439. 


±8. 


26.5 ± 0.45 


Leinertetal. (2001) 


1998 Jan 2 


1260 


±3 


-2.5 ±0.14 


432. 


±3. 


28.0 ± 0.45 


Leinertetal. (2001) 


1998 May 7 


1281 


±7 


-1.9 ±0.1 


408. 


±8. 


32.3 ± 0.73 


Leinertetal. (2001) 


1998 Oct 10 


1332 


±8 


-1.0 + 0.1 


377. 


± 19. 


38.2 ±0.5 


Leinertetal. (2001) 


1999 Jun 18 


1404 
1407 


±3 
±7 


0.44 ±0.1 
1.02 ±0.28 


318. 
303. 


±2. 
±7. 


49.4 ± 0.4 
52.1 ±0.3 


Leinertetal. (2001) 


1999 Aug 3 


Leinertetal. (2001) 


1999 Sep 1 


1414 


±4 


1.23 ±0.17 


292. 


±8. 


54.6 ±1.0 


Leinertetal. (2001) 


1999 Nov 23 


1437 


±4 


1.64 ±0.11 


279. 


±5. 


60.1 ±0.9 


Leinertetal. (2001) 


2000 Jun 4 


1487 


±9 


2.5 ±0.2 


240. 


±2. 


76.0 ± 0.4 


Leinertetal. (2001) 


2000 Jun 20 


1518 


±3 


3.6 ± 0.3 


237. 


±3. 


80.8 ± 0.7 


Leinertetal. (2001) 


2000 Oct 31- Nov 1 


1516 


±1.4 


4.4 ±0.1 


218.1 


±1.9 


88.9 ±0.6 


OCASS 


2001 Jan 9 


1546.4 


±1.1 


4.7 ±0.1 


219.0 


±9.0 


100.4 ± 2.3 


ADONIS, 66.C-0219 


2001 June 29 -July 6 


1600.7 


±5.5 


5.4 ±0.1 


222.7 


± 1.0 


120.3 ± 0.2 


SHARP, 67.C-0354 


2001 Dec 6 


1631.2 


±2.5 


6.8 ± 0.2 


244.9 


±0.4 


136.1 ±0.2 


NACO, 60A-9026 (commissioning) 


2001 Dec 11 


1633.0 


±7.7 


7.1 ±0.4 


246.7 


±3.8 


136.8 ±0.9 


ADONIS, 68.C-0539 


2002 Dec 16 


1712.1 


±3.4 


9.8 ±0.3 


327.7 


±0.9 


162.8 ±0.3 


NACO, 70.C-0476 


2003 Dec 12 


1772.3 


±4.9 


12.7 ±0.1 


406.3 


± 1.1 


178.2 ±0.1 


NACO, 072.C-0022 


2004 Dec 1 1 


1796.3 


±2.2 


15.9 ±0.1 


450.0 


±1.4 


189.7 ±0.2 


NACO, 074.C-0637 


2005 Sep 25 


1792. 


±19. 


17.9 ±0.2 


450. 


±11. 


198.7 ±0.9 


OCASS 


2006 Sep 1 - 2 


1759.2 


±4.6 


20.5 ± 0.2 


403.1 


±6.0 


208.1 ±0.9 


OCASS 


2006 Oct 30 


1742.9 


±2.9 


20.9 ± 0.5 


399.7 


±1.0 


209.6 ± 0.5 


NACO, 078.C-0386 


2006 Nov 28 


1744.7 


±3.8 


21.3 ±0.2 


396.1 


±2.0 


211.0 ±0.3 


OCASS 


2007 Sep 16 


1680.2 


±1.5 


23.6 ±0.2 


328.7 


±0.4 


223.5 ± 0.2 


NACO, 380.C-0179 


2008 Feb 1 


1646.4 


±1.4 


24.8 ± 0.2 


295.4 


±0.4 


231.5 ±0.2 


NACO, 380.C-0179 


2008 Oct 17,18 


1576.6 


±0.7 


27.3 ±0.1 


240.7 


±0.3 


252.6 ±0.13 


NACO, 382.C-0324/0329 


2008 Nov 6 


1571.3 


±0.7 


27.4 ±0.1 


237.6 


±0.3 


254.6 ±0.12 


NACO, 382.C-0324 



ing with 36.4 km/s to the Sun dBasri & Marcvill995h . Due to its 
kinematic properties, LHS 1070 can be associated with the old 
disk population. 

A first successful fit of the cl ose orbit of compone nt C around 
component B was presented by iLeinert et a l. (2001). Based on 
Speckle interferometric measurements a semi-major axis of a = 
446 + 29 mas and a period of P - 16.1 + L4yr were derived. 
When taking the refined distance of ICosta et al.l d2005l) into ac- 
count, the corresponding dynamical mass is 0. 157 + 0.042 M . A 
consis tent mass of 0. 157 + 0.009 M was found bv lSeifahrt et al.l 
(2008!!), who derived a wider (a - 461.9 + 0.7 mas) orbit, but 
with P - 17.0 yr. Due to the better coverage the authors pre- 
sented also a first reliable fit of the wide orbit of the components 
B and C around the primary. The dynamical mass of the whole 
system was determined as 0.272 ± 0.017 M . The orbits were 
found to be coplanar. 

In this study, we present new measurements taken both with 
Speckle interferometric techniques and adaptive optics in the K- 
band (Section 2). The data is carefully selected and calibrated on 
a case by case basis with the aim to reach the highest degree of 
consistency possible. These new data are not only used to fur- 



ther refine the fits of the close (Section 3) and the wide orbit 
(Sections 4). We utilise the primary component as astrometric 
reference, which enables the determination of the mass ratio of 
the components B and C. With this knowledge the individual 
masses of all three components can be derived. In Section 5 we 
discuss the implications of our findings with respect to the sta- 
bility of the system and its evolutionary stage. We conclude in 
Section 6. 



2. Observations and data reduction 

Our group monitored LHS 1070 since the discovery of its two 
companions, using a number of different telescopes and instru- 
ments. TableQ]gives a journal of observations and the measured 
relative positions. The following subsections describe the indi- 
vidual instruments and data reduction procedures. 



2.1. ESO NTT /SHARP 

The SHARP I camera (System for High Angular Resolution 
Pictures) of the Max-Planck-Institut for Extraterrestrial Physics 



2002 
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2004 2006 2008 



Fig. 1. Images of LHS 1070 obtained with NACO in December 2002, December 2004, October 2006, and November 2008. North is 
up, and East to the left. The separation between the two components B and C changes from about 330 mas in 2002 to about 240 mas 
in 2008 (cf. Table[Qi. 



dHofmann et all Il992h was already u sed fo r many of the ob- 
servations presented in iLeinert et al.l (1200 lb . We used it again 
during an observing campaign in June and July 2001 at the 
European Southern Observatory (ESO) 3.5 m New Technology 
Telescope (NTT) on La Silla, Chile. LHS 1070 was observed 
in the K-band at 2.2 /mi in four different nights. SHARP is a 
camera for speckle interferometry, which means that one ob- 
servation consists of several hundred frames with short integra- 
tion times (0.5 sec). To reduc e the data, we used our speckle 
program dKohler et al.. 2000). The same observing and data re- 
duction strategy was already employed for the data published in 
ILeinert et al. (2001 ) and a number of other multiplicity studies. 



2003). LHS 1070 was observed in the course of several pro- 
grams (see Table [TJ. To ensure a consistent data set, only 
imaging observations in the K s photometric band were used for 
the orbit determination. 

The NACO images were sky subtracted with a median sky 
image, and bad pixels were replaced by the median of the closest 
good neighbors. Finally, the images were visually inspected for 
any artifacts or residuals. Figure [T] sho ws an ex ample of the re- 
sults. The Starfinder program (Diolaiti et al., 2000) was used 
to measure the positions of the stars. The positions in several 
images taken during one observation were averaged, and their 
standard deviation used to estimate the errors. 



2.2. ESO 3.6m I ADONIS I SHARP II 

In January and December 2001, LHS 1070 was observed with 
the ad aptive optics (AO) system ADONIS ( Rousset & Beuzit, 
1 1999b and the SHARP II camera dHofmann et all Il995h at die 
ESO 3.6 m telescope on La Silla, Chile. We used the K-band 
filter of this instrument, which has a central wavelength of 
2. 177 /mi. The observing strategy was 'AO-assisted speckle in- 
terferometry", meaning we took many frames with short integra- 
tion times. In January, we recorded 1000 frames with an integra- 
tion time of 0.5 sec each. In December, we took 240 frames with 
a longer integration time of 3 sec each. For the data reduction, 
we used the same programs and algorithms as for the speckle 
data. 



2.3. CalarAlto 3.5m/ Omega Cass 

On several occasions between 2000 a nd 2006, LHS 1070 w as ob- 
served with the Omega Cass camera dLenzen et all 1 19981) at the 
3.5m-telescope on Calar Alto, Spain. The camera is equipped 
with a 1024x1 024-pixel detector, but we used only a subarray 
of 128x128 pixels to enable the fast read-out mode required for 
Speckle-interferometry. The observations were done in the K- 
band and with the highest resolution optics of the instrument 
(pixel scale ~95 mas/pixel). Our speckle program was used 
again to reduce the data. 

2.4. ESO VLT/ NACO 

The majority of new observations were taken with 
NAOS/CONICA (NACO for short), the adaptive optics, 
near-infrared camera at the ESO Very Large Te l escope (VLT) 
on Cerro Paranal, Chile dRousset et all 120031: iLenzen et all 



2.5. Plate scale and orientation 

For a heterogenous data set like this, it is crucial to calibrate the 
absolute pixel scale and orientation of each observation. To this 
end, we took images of fields in the Orion Trapezium during 
each observing campaign, and reduced them in the same way 
as the images of the science targets. The measured positions of 
the cluster stars were compared with the coordinates given in 
McCaugh rean & Staufferl (Il994l) . The mean pixel scale and ori- 
entation were computed from a global fit of all star positions. 
The scatter of values derived from star pairs were used to esti- 
mate the errors (see sect. 12.6) . 

The errors of the calibration are usually comparable to or 
larger than the errors of the measured positions of the science 
target, indicating the importance of a proper astrometric calibra- 
tion. For this reason, we decided to use in this work only data 
where images of the Orion Trapezium were taken within a few 
days. We have two observations of LHS 1070 that were taken in 
June and July, when Orion was not observable. These observa- 
tions are not used here, in order to ensure a consistent calibration 
of all the data. For the same reason, we do no use data from the 
literature that could not be re-calibrated by our group. 

2.6. Astrometric error estimates 

Independent of instrument and telescope used, our observation 
resulted in at least 4 images or speckle cubes for each epoch. 
The final result was obtained by averaging the relative positions 
measured in individual images. As estimate for the errors, we 
use the standard deviation. 

The error of the astrometric calibration was estimated in a 
similar way. We computed the pixel scale by dividing the sepa- 
ration of star pairs in arcseconds by the measured separation in 
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pixels. Similarly, the orientation is the difference of the position 
angle on the sky and on the detector. The standard deviation of 
the results for the many star pairs in our calibration field was 
used as estimate for the error. 

An additional source of error might be the motion of the tele- 
scope between the observation of the calibration field and the 
target. This should not cause much uncertainty for telescopes 
on a parallactic mount (i.e. the Calar Alto 3.5 m and the ESO 
3.6 m). However, telescopes on an alt-azimuth mount (ESO NTT 
and VLT) have to correct for image rotation, which might intro- 
duce an error on the order of 0.1° (W. Brandner, priv. comm.). 
Additionally, the opening and closing of the AO-loop, in con- 
nection with the active optics of the primary mirror, might in- 
duce small changes in the image scale. Our calibrations show 
that there are small changes in pixel scale and orientation on the 
time scale of weeks or months, but - to the best of our knowledge 
- there has been no study of short-term variations. 

This unknown error source might explain why we seem to 
underestimate our errors, which results in orbit fits with x 2 > 1 
(cf. tables0and|3). 



Table 2. Parameters of the best orbital solution for the pair B-C. 



Orbital Element 


Value 


Date of periastron To 


2454145.0 !j^ 




(2007 Feb 13) 


Period P (years) 


17 94+0-01 
1 ' - ZH -0.01 


Semi-major axis a (mas) 


457.8 33 


Semi-major axis a (AU) 


q o +0.07 
J.JJ _„ 07 


Eccentricity e 


fl f!997 +00004 
u.uzz./ _ 0000 3 


Argument of periastron a) (°) 


217 01 +004 
z-i /.ui _ 008 


P. A. of ascending node il (°) 


14.65 tg 


Inclination i (°) 


61-82 tgj 


System mass M s (mas 3 /year 2 ) 


(3.23 ± 0.01) • 10 5 


System mass M B + M c (M ) 


0.149 + 0.009 


Mean absolute difference 




observed - predicted (mas) 


4.9 


reduced^ 2 


3.5 



3. The orbit of LHS 1070 B-C 

We estimated the orbital parameters of the B-C pair by fitting 
orbit models to all observations listed in Table Q] We followed 
the procedure described in K ohler et al.l (I2008I) : a grid-search in 
eccentricity e, period P, and time of periastron Tq. At each grid 
point, the Thiele-Innes elements were determined by a linear fit 
to the observational data using Singular Value Decomposition. 
From the Thiele-Innes elements, the semimajor axis a, the angle 
between node and periastron u, the position angle of the line of 
nodes O, and the inclination i were computed. 

Since the orbit of LHS 1070 B-C is alrea dy quite well- 
known dLeinert et all 1200 it ISeifahrt et all 120081) . only a small 
range of parameter values had to be scanned: 200 points within 
0.018 < e < 0.028, 200 points within 17.21 yr < P < 17.31 yr, 
and initially 200 points for To distributed over one orbital pe- 
riod. After the initial scan over Tq, the best estimate for To was 
improved by re-scanning a narrower range in To centered on the 
minimum found in the coarser scan. This grid refinement was 
repeated until the step size was less than one day. 

We improved the results of the grid-search with a Levenberg- 
Marquardt y 2 minimization algorithm (Pr ess et aU [1992) that 
fits for all 7 parameters simultaneously. The simple-minded ap- 
proach would be to use the orbital elements with the mini- 
mum x 1 found with the grid-search. However, initial test runs 
showed that the algorithm does not converge on the global mini- 
mum. For the same reason, we did not use one of the previously 
published orb i t solu tions as starting point dLeinert et all 120011: 
Seifahrt et al., 2008). To make sure we find the globally mini- 
mum x 2 , we decided to use all orbits resulting from the grid- 
search as starting points that had x 2 < X 2 ■ + 9- The number 9 
was chosen arbitrarily, to avoid starting from obviously bad or- 
bits. The orbit with the globally minimum x 2 found in by the 
Levenberg-Marquardt fit is shown in Fig. |2] and its elements are 
listed in Table [2] To convert the semi-maj or axis from mas t o 
AU, we used the distance of 7.72 + 0.15pc (ICosta et all 120051) . 

The reduced x 2 of 3.5 is higher than expected for a good 
fit, which indicates that we underestimate our astrometric errors. 
The mean absolute difference between observed and predicted 
positions is also larger than typical errors of our measurements. 
This might be caused by unknown errors in the calibration, as 
mentioned in section l2~6l However, the residuals show no pattern 
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Fig. 2. The orbit of component C around component B. The 
observed positions are marked by their error ellipses and lines 
connecting the observed and calculated position at the time of 
the observations. The observations with NACO are marked by 
crosses. Their errors are too small to be discernible. The dash- 
dotted line indicates the line of nodes, the dashed line the peri- 
astron, and the arrow shows the direction of the orbital motion. 



indicating a systematic error, and the errors given in table [2] are 
our best guess. Therefore, we decided to accept the fit as it is. 

Errors of the orbital elements were determined by studying 
the x 2 function around its minimum. Since we are interested in 
the confidence interval for each parameter taken separately, we 
have to perturb one parameter (for example To) away from the 
minimum, and optimize all the other parameters. Any perturba- 
tion of a parameter will of course lead to a larger x 2 - The range 
in To within which x 2 (To) - X 2 ■ < 1 defines the 68 % con- 
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fidence interval for Tq. This interval is usually not symmetric 
around the To of the best fit, therefore we list in Table [2] sepa- 
rate limits for positive and negative perturbations. It should be 
noted that these limits describe the parameter range that contain 
68 % of the probability distribution, which is equivalent to the 
commonly used lcr-errors. However, the errors are not normally 
distributed, therefore a 2cr interval will not contain 95 % of the 
probability distribution. 

Estimating the error of the mass required a special pro- 
cedure. The mass itself is computed using Kepler's third law 
(M = a 3 IP 2 ). The semi-major axis a and the period P are usually 
strongly correlated. To obtain a realistic estimate for the mass er- 
ror, we did not use the naive way of error propagation. Instead, 
we considered a set of orbital elements where the semi-major 
axis was replaced by the mass. This is possible because Kepler's 
third law gives an unambiguous relation between the two sets of 
elements. With the mass being one of the orbital elements, we 
can treat it as one of the independent fit parameters and deter- 
mine its error with the method described in the previous para- 
graph. 

With an orbit derived from astrometric measurements, there 
always remains the ambiguity which of the two nodes is the 
ascending node (defined as the node whe re the companion is 
receding from the observer). Fortunately, Seifa hrtet al.l ([2008) 
measured the relative radial velocities of LHS 1070 B and C on 
2006 October 9. They found it to be negative, i.e. LHS 1070 C 
was approaching u£]. Therefore, the position angle of the as- 
cending node is 14.65° (cf. Fig.|2]i, and the position angle of the 
descending node is 194.65°. Our orbit model predicts a relative 
radial velocity of -4.8 ±0.1 km/s on 2006 Oct 9, in agreement 
with the measurement of Sei fahrt et al.l (2008). 



4. The orbit of LHS 1070 BC around A 

The orbit of components B and C around each other allows us to 
determine only the combined mass of B and C. To compute the 
individual masses, we need to know the mass ratio q, which can 
be computed if the position of the center of mass (CM) of B and 
C is known. Unfortunately, we cannot observe the CM directly. 
However, we know that the CM of B and C is in orbit around 
component Aq and that B and C are in orbit around their CM. 
The CM is always on the line between B and C, and its distance 
from B is the constant fraction q/(l + q) of the separation of B 
and C. 

We follow the method that was used bv lKohleret alJ J2008I) 
to derive masses in the triple system T Tauri. The position of 
the CM of B and C is described in two ways: First, it is on a 
Keplerian orbit around A, which is described by 7 orbital ele- 
ments. Second, the position of the CM can be computed from 
the observed positions of B and C, and the mass ratio (which is 
treated as a free parameter). Standard error propagation is used 
to obtain an error estimate for this position. To compute x 2 , we 
compare the position of the CM from the orbit around A with the 
positions derived from the observations. Our model has therefore 
8 free parameters, the 7 elements which describe the orbit of the 
CM of B+C around A, and the parameter / = ql{\ + q ). The 
parameter / is often called fractional mass (Heintz] .[l978l) . since 



Table 3. Parameters of the best solution for the orbit of BC 
around A. 



Orbital Element 






Value 


Date of periastron To 






2459727.9 *** 
(2022 May 28) 


Period P (years) 






44 4+11-9 
^*" -2.4 


Semi-major axis a (mas) 






1111.6!™ 


Semi-major axis a (AU) 






o co +0.17 
8.D8_ 017 


Eccentricity e 






5200 +0002 ° 

U - Jiuu -0.0012 


Argument of periastron w 


(°) 




147.61 ?™ 


PA. of ascending node Q 


(°) 




26.80 !^ 


Inclination i (°) 






54 75 +Q-90 

J J -0.10 


System mass M A + M B + 


M c 


(mas 3 /yr 2 ) 


(6.98 ± 0.02) ■ 10 5 


System mass M A + M B + 


M c 


(M e ) 


n a 7 i +0.020 


Mass ratio M C /M B 






0.92 ± 0.01 


Mean absolute difference 








observed - predicted (mas) 


7.9 


reduced x 1 






2.6 


Mass of A M A (M Q ) 






0.172 + 0.010 


MassofBM B (M ) 






0.077 ± 0.005 


Mass of C M c (M Q ) 






0.071 ±0.004 



1 It is not entirely clear whether Seifahr tet al.l (12008b give the veloc- 
ity of C minus B, but they report £2 = 14.5°, which is consistent with 
our interpretation 

2 We did not detect the companion D reported bv lHenrv et al.llll999h . 
Therefore, by "A" we mean the suspecte d close binary com posed of A 
and D (if it exists) in the nomenclature of lHenrv e t al. ( 1999). 



it is the secondary star's fraction of the total mass in a binary. It 
is useful in our case because it also describes the fractional off- 
set of the CM from B, i.e. the separation between B and the CM 
divided by the separation between B and C. For a grid-search, / 
is better suited than q, because / is confined to the range to 1, 
while q is a number between and infinity. 

The fitting procedure is similar to that used for the orbit 
of B-C, except that the grid-search is carried out in 4 dimen- 
sions: eccentricity e, period P, time of periastron To, and the 
fractional mass /. Singular Value Decomposition was used to 
fit the Thiele-Innes constants, which give the remaining orbital 
elements. It is worth noting that the orbital elements in this fit 
describe the orbit of the A-BC binary, only the fractional mass / 
refers to the pair B+C. 

For the four dimensional grid search, we used 100 points 
each for /, e, and P. The grid ranged from to 0.99 in /, to 
0.99 in e, and 30 to 300 years in P. The grid in Tq started with 
100 points distributed uniformly over one orbital period. Similar 
to the fit for the orbit of B-C, the grid in 7o was refined until the 
grid spacing was less than one day. 

The results of the grid-search were impro ved with a 
Levenberg-Marquardt x 2 minimization algorithm dPress et all 
1992). As starting points, we used all orbits resulting from the 
grid-search that had x 2 < X 2 ■ +9. The orbit with the globally 
minimum x 2 found by the Levenberg-Marquardt fit is shown in 
Figs.|3]and|4] Tablef3]lists its elements, as well as the mass ratio 
q, and the individual masses of all three components. As with 
the inner orbit of B and C, x 2 is somewhat larger then expected 
for a good fit. 

Errors for the parameters were again estimated by analyz- 
ing the x 2 function around the minimum. Since only a relatively 
small fraction of the orbit has been observed so far, the uncer- 
tainties for the orbital elements are much larger than for the or- 
bit of B and C around each other. Figure |5] shows x 2 as func- 
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Fig. 3. The orbit of the center of mass of components B and 
C around component A. The observed positions are marked by 
their error ellipses and lines connecting the observed and cal- 
culated position at the time of the observations. Note that these 
positions are computed from the observed positions of compo- 
nents B and C, and the mass ratio q, which is a free parameter 
of the model fit. The dash-dotted line indicates the line of nodes, 
the dashed line the periastron, and the arrow shows the direction 
of the orbital motion. The two orbits shown are the model with 
the globally minimum x 1 , and an almost circular orbit with a 
longer period of about 80 years. 



tion of period and eccentricity. Orbits with periods of more than 
100 years are within the 95 % confidence region, and periods of 
more than 200 years are in the 99.7 % confidence region. The er- 
ror of the system mass was derived in the same way as with the 
orbit of LHS 1070 B and C around each other, i.e. by treating M 
as an independent fit parameter. 

The ambiguity of the ascending node is resolved with the 
help of Seifahrt et al. (2008). They measured in 2006 a positive 
relative radial velocity between the barycentre of LHS 1070 B 
and C and component A. Therefore, the binary BC was receding 
from the observer, and the position angle of the ascending node 
must be 26.8° (cf. Fig. |3). Our orbit model predicts a relative 
radial velocity of 3.1 + 0.1 km/s on 2 006 Oct 9, in agreement 
with the measurement of Seifahrt et al.l (120081) . 

Finally, the mass of component A can be computed from the 
system mass of the triple and the mass of the binary B+C (de- 
rived in section [3). The largest contribution to the error of both 
masses is the uncertainty in the distance to the system. Since 
this contribution to the errors of Mabc and Mbc is correlated, it 
would not be correct to add the errors in quadrature. Instead, we 
computed masses from the semi-major axes in mas (which re- 
sults in the unusual mass unit of mas 3 /years 2 ), subtracted Mbc 
from Mabc, and converted the resulting Ma into solar masses by 
multiplying it with the distance cubed. This way, the error of the 
distance enters the calculation only once, resulting in the correct 
estimate for the error of Ma- 



Fig. 4. Enlarged section of the orbits in Fig. [3] showing the part 
covered by observations. 




Fig. 5. x 1 as function of P and e for the orbit of BC around A. 
The cross at P — 44.4 yr, e - 0.52 marks the minimum. The 
contour lines encircle the 68 % confidence region (which is a se- 
ries of mostly unconnected patches), the 95 %, and the 99.7 % 
confidence region (corresponding to lcr, 2cr, and 3<x in the case 
of normally distributed errors). The dashed line marks the stabil- 
ity limit for the system, only orbits below this line are long-term 
stable (see section IBTTT i. 



5. Discussion 

5.1. Is the LHS 1070 system stable? 

Our best fit for the orbit of the binary BC around A is only a 
factor 2.5 larger than the orbit of BC itself. The outer orbit has 
also a rather large eccentricity. As a result, the distance from B 
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or C to A becomes smaller than the distance to its binary com- 
panion during the periastron passage of the outer orbit (Fig.[6]). 
This raises the question whether a triple system like this can be 
stable over timescales comparable to its age. 




2000 

Year 



Fig. 6. Distances between the three stars in the system as func- 
tion of time. Here we assume that the stars stay on the Keplerian 
orbits derived in sections [3] and |4] i.e. we do not simulate the 
three-body problem posed by the triple system. 
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Donnison & Mikulskis (1995) presented a criterion for the 
long-term stability of coplanar triple systems: The periastron dis- 
tance qi - «2(1 - £2) °f tne oute r orbit should be at least a factor 
of 3.3 larger than the semi-major axis ci\ of the inner orbit. In 
non-coplanar systems, increasing the inclination decreases the 
region of stability , i.e. th e outer orbit has to be even larger to be 
stable (Donnison, 2009). In Fig. [5J the limit for coplanar orbits 
(q2 > 3.3 ai) is indicated by the dashed line, only orbits below 
the line are stable. The orbit with the minimum^- 2 is clearly in 
the unstable regime. However, the 2<x and 3cr confidence regions 
extend well into the regime of stable orbits. It is therefore possi- 
ble that the true orbit has a longer period of at least 80 years (cf. 
Fig.E). 

How does this uncertainty about the orbit influence our con- 
clusions? Figure |7] shows the parameters of the outer orbit that 
we are most interested in. Shown are all orbits within the 99 % 
confidence region in our grid of possib le solutions. Unstable or- 
bits (according to the criterion bv lDonnison & Mikulskis, 1995) 
are indicated by gray symbols, stable orbits are black. The top 
panel shows the angle between the planes of the inner and 
outer orbits. This angle is 12.5° for the orbit with the mini- 
mum^ 2 , and less than 10° for all stable orbits. We find a few 
orbital solutions with angles larger than 100°, but all these or- 
bits are unstable. We therefore conclude that the inner and outer 
orbit are almost coplanar. Th erefore, the stability criterion by 
Don nison & Mikulskis! (1 19951) is sufficient in our case, although 
it does not apply to inclined orbits. 

The second panel of Fig. [7] shows the masses of LHS 1070A 
and B. Orbits with short periods result in unrealistically high 
masses for A, which also indicates that the true period is longer. 



Fig. 7. Results of the fits for the outer orbit as function of 
their period. Shown are all model orbits within the 99 % con- 
fid ence region (3<x). Unstable or bits (according to the criterion 
bvlDonnison & Mi kulskis! 1 19951) are indicated by gray symbols, 
stable orbits are black. The orbit with the minimal x 1 is marked 
by an asterisk. 



The mass of component B lies between 0.076 and 0.077 M , 
almost independent of the orbital period. The mass of com- 
ponent C was not plotted, since it can easily be computed as 
O.149M -M B . 

Finally, the bottom panel of Fig. |7]shows the mass ratio be- 
tween LHS 1070 C and B. All stable orbits result in a mass ratio 
of 0.923, which is the same as the mass ratio of the orbit with 
the minimum^ 2 . 

We conclude that it is possible to find stable orbital solutions 
that are compatible with the astrometric measurements. In these 
stable configurations, the orbit of B and C around each other, and 
the orbit of B+C around A are almost coplanar, with angles of 
less than 10°. The mass of component A lies between 0.13 and 
O.16M for stable orbits. 

As a final note on the stability of the system, we point out 
that the stability criterion used here does not take into account 
resonances. For example, a period of the outer orbit of 52 years 
would be in a 3: 1 resonance with the inner orbit. This orbit would 
be within the \cr region around the orbit model with minimum 
X 2 - However, the uncertainties of our orbit fit are too large to pro- 
vide a reasonable starting point for a search for stable resonant 
configurations. Furthermore, it would not significantly change 
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our conclusions about coplanarity and the masses of the three 
components if the period of the outer orbit was between about 
50 and 80 years. 

5.2. Comparison with theoretical models 

The goal of dynamical mass determinations is to test theo- 
retical predictions for the mass and luminosity of the stars. 
Figure H] shows a mass -luminosity-diagram with the isochrones 
of iBaraffe et al.l d 19981) and our results for the compone nts o f 
LHS 1070 (using the V-magnitudes given by Leinert et al. 2000). 
Ad ditionally, it shows the empirical Mass-Luminosity-Relation 
by iHenrv et all d!999l) . which is similar to the theoretical 500- 
Myr-isochrone. Within the error bars, both LHS 1070 B and C 
are compatible with an age between 500 and 800 Myr. 

On the other hand, LHS 1070 A is either 3 mag too faint or 
0.07 M too heavy compared to the theoretical models and the 
empirical relation. This is also true if we take into account that 
stable orbit configurations result in lower masses (indicated by 
the dotted line in Fig.©. In the K-band, LHS 1070 A is about 
l mag too faint compared to the model isochrones. This excludes 
the possibility that the under-luminosity is due to extinction by 
circumstellar or foreground material (as unlikely as it may be 
for a star of this age and distance), since then the effect in the 
infrared should be smaller, only about 0.3 mag . 
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Fig. 8. Mass-luminosity diagram with the 
LHS 1070, the models of Baraffe et al 



components of 
(119981). and th e empir- 
ical Mass-Luminosity-Relation by [Henry et all (Il999l) . Shown 
are the theoretical isochrones for ages between 500 Myr and 
5 Gyr. The dotted horizontal line indicates the mass range for 
A resulting from stable orbit configurations. The points labeled 
'Aa photometry" and "D photom etry" indicate the m agnitudes 
if A is split into two components (IHenrv et all [T.999), with the 
masses adjusted to put the points onto the 630-Myr isochrone. 



estimate the mass of Aa and D to be 0.10M o and 0.08 M , resp. 
The combined mass of A resulting from photometry is there- 
fore about 0.18 M , in reasonable agreement with the dynamical 
mass resulting from the orbit determinations. 

In an attempt to find a possible close companion to 
LHS 1070 A, we analyzed all available high-resolution data of 
the system. The resolution required to detect component D can 
only be reached with a 8 m-class telescope, preferentially at short 
wavelengths, e.g. J-band. This means we have to rely on NACO 
data. There are two datasets taken with NACO in the J-band 
on June 27. and December 12., 2003, and two datasets in the 
narrow-band filter NB_2.17 in a mode that allows speckle post- 
processing of the data. In all four datasets, we can exclude an 
companion with a flux ratio brighter than 0.2, i.e. a magnitude 
difference of 1.75 mag . The resolution of the four observations is 
30 - 40 mas. With the masses of Aa and D derived above, the 
magnitude difference should be about 0.8 mag . Assuming a semi- 
major axis of 50 mas and a system mass of 0.17 M results in 
an orbital period of about 0.6 years. It is possible that all our ob- 
servations were carried out at times when the companion was 
too close to the primary to be visible. However, given the large 
number of attempt to detect LHS 1070 D this appears unlikely. 

We also looked into the possibility to detect a periodic astro- 
metric signal of LHS 1 070 D in the residuals of the A-BC separa- 
tion. If LHS 1070 A is an unresolved binary, then an astrometric 
shift of its center of light with the period of its orbit can occur. 
The magnitude of this shift depends on the position of the center 
of mass and the position of the center of light, hence the mass 
ratio and the flux ratio. The mass ratio of our suspected Aa-D 
binary is 0.8. while the flux ratio computed from the models of 
IBaraffe et all (Q998) is 0.48. The expected motion of the cen- 
ter of light is therefore about 0.12 times the size of the orbit of 
D around Aa (6 mas for an orbit of 50 mas). Unfortunately, our 
sampling is not very dense, with 38 observations in 15.25 years. 
The Nyquist frequency v^ = N/(2T) corresponds to a period 
of 0.8 years. Since our sampling is rather uneven, the Nyquist 
frequency is not a sharp limit for the detectability of a periodic 
signal. Nonetheless, due to the sparse sampling, it is unlikely 
that we can detect periods much shorter than 0.8 years. In fact, 
neither classical nor generalized Lomb-Scar gle-Periodograms 
dScargld Il982t IZechmeister & Kiirsterl 120091) show any signif- 
icant peaks. Even adding an artificial signal with a period of 
0.6 years and a large amplitude of 100 mas does not result in 
a significant detection, although peaks at 1/2 and 1/4 of the fre- 
quency appear. In summary, we conclude that an astrometric de- 
tection of LHS 1070 D is not possible with the available data. 

Another possible explanation could be that LHS 1070 A is 
a binary too close to be resolved by the observations available 
so far. This hypothetical companion would not be identical with 
component D. However, at the moment the most likely explana- 
tion might be that we overestimate the system mass. This would 
not be too surprising, given the small observational coverage of 
the orbit. 



The simplest explanation would be that A is indeed a bi- 
nary as reported by IHenrv et all (1999). In the following, we 
will call the components Aa and D, and continue to use A for 
the (unresolved) binary system. We computed the absolute V- 
magnitudes of Aa and D from the c ombined apparent magnitude 
V = 15.35 mag (Leinert et al., 2000) and the magnitude difference 
AV - 2.46™ g (Henry et al., 1999), using the distance of 7.72 pc 
dCosta et all 120051) . By folding th e intersection with the 630- 
Myr isochrone from the models of Baraffe et al. ( 1998), we can 



6. Summary and Conclusions 

We present new relative positions of LHS 1070 A, B, and C, col- 
lected in the years 2000 to 2008. They were used to derive an 
improved model for the orbit of B and C around each other, 
and an estimate for the orbit of B and C around A. The orbit 
of B and C is well-determined by now, with an orbital period of 
17.24 + 0.01 years and a system mass of 0.149 + 0.009 M (cf. 
Table©. 
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The observations of LHS 1070B relative to A span only a 
range in position angle of 36°. The orbit of B and C around 
A derived from these observations is therefore not so well- 
constrained. The orbital solution minimizing x 2 nas a rather 
short period of about 44 years, but it results in an unstable con- 
figuration of the triple system. It is more likely that the true or- 
bital period is in the range 80 - 200 years. Orbital solutions in 
this range of periods are still within the 95 % and 99.7 % con- 
fidence region for the fit to the astrometric observations (which 
corresponds to 2<x and 3<x for the case of normally distributed 
measurement errors). Despite this uncertainly, we can constrain 
the system mass M A + M B + M c to 0.28 to 0.31 M if we accept 
only stable orbit configurations. 

The outer orbit also yields the mass ratio of B and C, which 
is quite well-constrained to 0.92 + 0.01. Taken together, these 
results yield individual masses of Ma = 0.13 . . .0.16M o , Mb - 
0.077 + 0.005 M and M c = 0.07 1 + 0.004 M , 

Placing the three stars in a mass-luminosity diagram and 
comparing with theoretical isochrones shows that B and C are 
coeval within the measurement errors, with an age between 500 
and 800 Myr. On the other hand, LHS 1070 A appears to be too 
faint for its mass, or too massive for its luminosity. One pos- 
sible explanation could be that LHS 1070 A is itself a binary. 
There has been one report of the discovery of a close companion 
dHenrv et all 11999) . but this discovery could not be confirmed, 
despite the large number of observations collected. The detec- 
tion might have been caused by a glitch in the data (T. Henry, 
priv. comm.). To confirm or disprove the binarity of LHS 1070 A, 
observations sensitive to companions at very small separations 
would be helpful, either spectroscopic or interferometric . For the 
time being, the most likely explanation might be that our orbit 
fit overestimates the system mass. 

Our results for the masses of the three stars show that 
LHS 1070C is almost certainly a brown dwarf. LHS 1070B is 
very close to the hydrogen-burning mass-limit, possibly also a 
brown dwarf. While LHS 1070 A is clearly above the hydrogen- 
burning limit, it is still a very low mass star. The very low 
mass triple system LHS 1070 thus has the very interesting prop- 
erty to contain one component above, one below, and one 
at the hydrogen burning limit. This combination makes it an 
ideal candidate for testing the change of atmosphere proper- 
ties in this mass regime. Furthermore, it allows some conclu- 
sions about its formation. One theory for the formation of brown 
dwarfs is that they got ejected out of a multiple system be- 
fore tiiej^_could_accrete enough mass to start hydrogen burn- 
ing (JReipurth & Clarkel 1200 ll) . It is difficult how a system like 
LHS 1070 can survive such an event. Even if it managed to re- 
main bound, one would expect more eccentric and inclined or- 
bits than found in this work. It is more likely that LHS 1070 
formed in the same way as normal stars, where a number of triple 
and higher-order multiple systems are known. 
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